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Abstract 

We study multivariate generalizations of the g-central limit theorem, a generalization of the 
classical central limit theorem consistent with nonextensive statistical mechanics. Two types of 
generalizations are addressed, more precisely the direct and sequential (/-central limit theorems 
are proved. Their relevance to the asymptotic scale invariance of some specially correlated 
systems is studied. A g-analog of the classic weak convergence is introduced and its equivalence 
to the q- convergence is proved for q > 1. 

1 INTRODUCTION 

We studied (/-analogs of the classical central limit theorem and a-stable distributions for ^-independent 
random variables in [TJ [H [3] . The results that we obtained there are consistent with nonextensive 
statistical mechanics [U E] and represent limit states of sequences of globally correlated systems. 
These limit states are described by g_i-Gaussian distributions, where 

3a- 1 

9-1 = — rr > 
<7 + l 

where the index q characterizes the correlation. The particular case q = 1, for which g_i = 1, 
corresponds to independent random variables, and exibits the classic central limit theorem. If 
g / 1 then g_i / 5 and there exist series of g-CLTs depending on the type of correlation (see 
PQ). The main tool used for the study of the g-central limit theorems is the q- Fourier transform 
introduced as 

F q [f] = J f{x) ® q e^dx = J fWeftVW-'dx, 

where / : R — > R\ is a density function, and e q {x) is the ^-exponential (for elements of the (/-algebra 
used in the current paper, see [9]). F q is a nonlinear transfomation if q / 1. 

In the present paper we will study a multivariate analog of the g-CLT. We will introduce two 
generalizations, denoted direct and sequential, of the (/-Fourier transform to the multidimensional 
case. In turn, they lead to the two generalizations of the g-CLT that we address here. Finally, we 
apply the multivariate g-CLT to the study of the scale invariance of statistical quantities with a 
finite g-mean and arbitrary (2q — l)-variance matrix. 
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2 PRELIMINARIES 

Let (Q,J-,P) be a probability space and X : Q — > i? d be a d-dimensional random vector with a 
joint density function fx( x ) = /(^)> x = (^li •••> %d) £ We introduce the associated escort 
density |10j 

M ' ~ v q {f) ' 

where v q {f) = [f{x)] q dx. By definition, the (/-mean of X is a vector Hq(X) = fj, q = (nq t x, ^ q ,d) T , 
where = fj, q j(Xj) = J Rd Xjf q (x) dx, j = 1, .., d, and the (/-covariance matrix S g is a matrix with 
entries = Spd{xj — fi q j)(xk — fJ- q ,k)fq(x)dx. By the direct d-dimensional (/-Fourier transform 
we understand the integral 

F q [fKO= [ m® q e%* dx= j fixyf^-'dx, (1) 

JR d JR d 

where £ £ R d and ix^[f(x)] q ~ 1 = i[f(x)] q ~ 1 J2k=i x kiki subject to the existence of the integral. Let 
q = 1. Then, (JTJ) identifies the classic Fourier transform 



m(0 = I f(x)e^dx 

JR d 



R d 

Moreover, in this case J- = F\ o ... o F d , where Fj, j = l,...,d, are the one-dimensional Fourier 
transforms with respect to the variables Xj, j = 1, d, and the symbol "o" means the composition 
of operators Fj, j = 1, ...,n. This property fails if q ^ 1. That is, generically 

Fq / Fl,q ° ••• ° F dtq , 

where Fj >q is the (/-Fourier transform with respect to Xj, j = 1, ...,d. Assume 

_ 2g + n(l - q) 
qn ~ 2 + n(l-g) 

i.e., equivalently, 



n = 0,±l,±2,... (2) 



1 1 n 

+ -, n = 0,±l,±2,... (3) 



1- q n 1-q 2 

Let Qk be an ordered set of numbers Qk = {qk, Qk+li •••) Qk+d-l}- The composition 

T Qk [f](0 = F ltqk+d _ 1 o...oF dtqk [f(x 1 ,...,x d Mu--^d) (4) 

is called a Qk- sequential Fourier transform of a given function f(x). The direct and sequential 
Fourier transforms coincide, only if q = 1. In this case Q = {1, 1}. 
A d-dimensional (/-Gaussian with a covariance matrix £ is a function 

G q (X;x) = C d , q (P)e q ^ x \ (5) 

, positive real number, 

1, d, and 



where (3 is a positive real number, XI is a positively defined d x (i-matrix with entries — , i,j 



C d , q {P) = C d! q(f3, S) = ( / e-f&rtdx 

R d 



is the normalising constant. We denote C<j )9 = C<j )? (l,£). Clearly, C djq ((3) = C d ^ q {l, /?£). In the 
particular case of a diagonal covariance matrix £ with diagonal entries Oj.j > 0, the (/-Gaussian 
takes the form 

G q (£;x) = C d , g (/3)e/ (ffl ' l|xi|2+ "- +ffd ' d|Xd|2 \ 
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Now we compute the value of C d ^ q {(3) as a function of q, (3 and X. First, assume (3 = 1 and 
dei(E) = |S| > 0. There is an orthogonal matrix U such that £ = U T DU where D is a diagonal 
matrix. Substituting y = Ux in the integral {C d-A )~ 1 = J^d e g ^'^dx we have 

(C dq y l =( e -^dy = -±=( e-^dz = ^4, 

d 

where oj d = ^nk is the area of the surface of the unit sphere in R d and 



roc „ 

I q ,d = / r*~V dr. (6) 



o 



Note that <ooifq , <l + §. Hence, 



= (7) 



We also note that consequently integrating the d-dimentional g-Gaussian, for C d-q with T, = I 
(I is the identity matrix), we obtain the formula 

(3_^(3_^...(: 



Lyg 9l ...U gd _ 1 

with q, qi, q^-i identified in ([2]), and C q , the normalising constant of the one-dimensional Gaus- 
sian. 

For an arbitrary (3 > 0, ([7]) implies the following lemma. 
Lemma 2.1 For i/ie normalising constant of a q-Gaussian = C dt qe q ^^ x ' x \ the formula 

C d M = PiV f m =^C d ,q (9) 
^ 2 Iq,d 

holds, where T(-) is the Euler's gamma-function and I qd is defined in 

Lemma 2.2 Let q < 1 + 1 and S 6e a symmetric positively defined dxd matrix with a determinant 
|E| > 0. T/ien 

2-d(g-l) 

' . v I , , , \ 

d,q 



e q ^™ I . (10) 



V 



Proof. First we assume that (3=1. Applying elementary properties of ^-exponential and 
g-product we have 

UC d , q e~ { ^ x) m = C d , q I eq {Sx ' X)+lX ^ dx. (11) 

JR d 

Further, by assumption, £ is a symmetric positive matrix. Therefore there exists an orthogonal 
matrix U such that D = [7 _1 EC/ is a diagonal matrix. Assume S±, ...,S n are the diagonal entries of 
D. Then (fTTj) is reduced to 

r,lcW^U) = %jf «f «. ^' {D -' U ^, (12) 
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where |S| = 5\...8 n is the determinant of X. Again using the properties of ^-exponential and 
g-product we compute the right hand side of (fl~2j) : 



Finally, for 0^1, replacing £ by (3T, and taking into account ([9]), we arrive at (jlOp . 
Corollary 2.3 Under the conditions of Lemma \2.2\ 

F q lC d ,Me« mx ' x) }(0 = e- a{ ^\ (13) 

w/iere 

_ 2g-d(g-l) _ 2(g-l) 
91 " 2-d(g-l) - 1+ 2-d(g-l) (14) 

and 

(2 - rf(g - l))[C d Mf^) 

8/3 1 J 

Replacing Cd, q {(3) by its value in equation Q we have 



a 



2-d(g-l) /^(j)V M> 



2 2g+l / gl+d-d 97r d(9-l) I J d 
Introduce the function 



2s-d(s-l) n , 2(a-l) 



A«) = o ,/„ TV = 1 + 



2-d(s-l) 2-d(«-l)' 
and denote its inverse which reads 



, .._2t + rf(t-l) , , 2(i-l) 



= - 1 r^r — ^ = i + 



2 + d(t-l) 2 + rf(t-l) 

z(s) possesses the property 



One can easily verify that z d (j^^) = £ and z d (^) = z d,-i^ ■ Moreover, it is readely verified that 



(z^s))- 1 +dz d '- 1 (s) = l + d. (16) 

Further, let us introduce the sequence q d = z d (q d _ 1 ),n = 1,2,..., with a given go = g, g < 
1 + ^ ■ We can extend the sequence q n for negative integers n = —1, —2, ... as well, putting g_ n = 
z~ (qf- n ), n = 1, 2, ... . It is not hard to verify that 

d _ 2g + rfn(l-g) 2(g-l) _ 

9 » 2 + dn(l-g) + 2-dn(g-l)' " " U ' ±L ' ±2 ' ~ U0 

i.e., equivalently, 

1 1 dn , . 

+ — , n = 0,±l,±2,... (18) 



1 - ^ 1-9 2 

Note g^ = 1 for all n = 0, ±1, ±2, if g = 1 and lim n ^ r ± 00 q li l = 1 for all g ^ 1. 
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Remark 2.4 Essentially the same mathematical structure has already appeared in a quite different, 
though possibly related, context: see Footnote of page 15378 of In the one- dimensional case, 
i.e for d = 1, the relationship (T7\ ) is obtained in JT$. For n = 1, in the d-simensional case, the 
relationship 

f= 2g-< g -l) 2( g -l) 

between components of elliptically invariant Gaussians was recorded in Ji£J /. 
Lemma 2.5 The members of the sequence qf t satisfy the following duality relations 

d^_i + -^- = l + d,n = 0,±l,±2,.... (20) 

Proof. By elementary calculations we have 

n , , d 1 2d(l -q) 1 1 

1 + d- dql_ x = 1 + - 



2 + d(n-l)(l-g) 2 g +d(n+i)(i- g ) „d ■ 

V A y; 2+d(n+l)(l-g) n+1 



d(n+l)(l-g) 

Remark 2.6 In the one- dimensional case (d = 1) the duality relations A20\) recover the well known 
relationships (see, e.g., fl\ [13, \1J$) 

q n -i H — = 2, n = 0, ±1, .... 

?n+l 

It follows from Corollary 12.31 and the definition of the sequence q% the following statement. 
Lemma 2.7 The series of mappings 

F q d F q d F q d j. F q d F q d 

... Q q d_ 2 Q qd _^ -3- G q -4 g qd -4 g qd -4 ... (21) 

^ jC-1 %d %d 

... ^ 3 a g d 2 ^ 2 <- a 9 a g g d g qd ... (22) 



hold. 

Lemma 2.8 



F q [C Q e-^+ A \Z) = Cxe-A^W"), (23) 

■9-1 



twffc Ci = /3i = and a x = ^ 

Corollary 2.9 Ze£ /3j > 0, j = 1, ...d, 6e given numbers. Then 

„9fe+d-l -1 
_y^d Td-j + l C d _j 2 

J~Q k [(-'Qe qk \ - C d e qk+d , (24) 



3=1 j=m+l 

and 



Remark 2.10 It follows from |2j) that, if we take Co = ^^"'^ d fin which case Coeql Xl+ "' + ^ dXd is 
a Gaussian density), then Cd = 1 in (24\ ). 
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3 THE g-CONVERGENCE 



In this section we introduce the notion of the weak q-convergence of a given sequence of random 
variables, and compare it with the multivariate generalization of the q-convergence. A sequence 
of random vectors is said to be g-convergent to a random vector X if .F g [Xjv](£) — > ^[X]^) 
for every £ G K, where K is arbitrary compact in R d . In the one-dimensional case, this definition 
was introduced in [I]. In terms of associated densities the g-convergence is equivalent to the locally 
uniform (by £) convergence ^ [/tv] (<!; ) — ► [/](£)> where fjy and / are densities of Xjy and X 
respectively. We denote the q-convergence by the symbol — >. 

Recall that a sequence Xn is called weakly convergent to a random vector X if the sequence of 
corresponding distribution functions Fx N (x) converges to Fx(x) at every point of continuity [15]. 
The weak convergence is usually denoted by =>. The equivalent definition of the weak convergence 
in terms of associated densities is the following: for any bounded continuous function g(x), the 
sequence of integrals J R fN(x)g(x)dx converges to f R f{x)g{x)dx |15| . For q / 1 denote by W q 

the set of functions <f> continuous on R d and satisfying the condition \4>(x)\ < C(l + x G 

R d . A sequence of random vectors Xjy is called weakly (/-convergent to a random vector X if 
J Rd fN(x)dm q — ► J Rd f{x)dm q for arbitrary measure m q defined as dm q () = (p q (x)dx, where (p q G 
W q . We denote the weak ^-convergence by the symbol Notice that, in the 1-dimensional case 
{d = 1), any measure m q , 4> q G W q is finite (i.e. m q (R 1 ) < oo) for all q > 1. If d > 2 then for 
1 < q < 1 + l/((i — 1) a measure m g is finite (i.e. m q (R d ) < oo). If < q < 1, then there exists a 
function G W g such that the corresponding measure m q diverges (i.e. m q (R d ) = oo). 

Lemma 3.1 If q > 1 then Xn =$> Xq yields Xn -^o- If Q < 1 then Xjy =l> yields Xn Xq. 

The proof of this lemma follows from the obvious fact that W q is a subset of the set of bounded 
functions if q > 1, and W g contains the set of bounded functions as a subset if q < 1. 

The next lemma uses the classic notion of tightness of a sequence of probability measures. 
Recall that a sequence of probability measures is called tight if, for an arbitrary e > 0, there is 
a compact if e and an integer N* such that HN(R d \ K e ) < e fox all iV > N*. 

Lemma 3.2 Lei 1 < q < 1 + l/(d — 1) if d > 2 and l<q<2ifd = l. Assume a sequence of 
random vectors Xn , defined on a probability space with a probability measure P, is q-convergent. 
Then the sequence of associated probability measures = 

P{X~ 1 ) is tight. 

Proof. It suffices to prove the statement in the one-dimensional case. The multivariate extension 
can be obtained using the analogous arguments proceeded in the classic case (i.e. in the case q = 1) 
using the Cramer- Wold device (see, for instance, [HI US]). Thus, assume that 1 < q < 2 and 
Xjy is a ^-convergent sequence of random variables with associated densities /at and associated 
probability measures fiN- We have 



1 



/ (1-F ff [/*](0R = ~ / (1- / fNef^ 1 dx)di 

J-R R J-R JR 



R 




(26) 



It is not hard to verify that 




(27) 
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It follows from ([26]) and ([27]) that 

i r R roc ( sin^_(x{2 - q)Rf^~ 1 ) . 

r I J 1 - F ' lf " m)di = 2 L 1 — fcB-rirtr 1 <28) 



iV 

Since 1 < q < 2 by assumption, > 1 as well. It is known [2] that for any q > 1 the properties 
sin j(x) < 1 and (sin '(x))/x — > 1, x — > hold. It follows from these facts and {x : (a;) > 

fl(2-g) } C i X 1 N > R(2-q) } tllat 

Further taking into account the g-convergence of /jv to /o and taking i? small enough we obtain 

-R 
-R 



UN (\x\ >r)<^ J\l ~ F q [f Q }(0)<% < e, 



where r = ji^-q) » anc ^ 6 * s an arbitrary positive number. 
Introduce for t G [0, 1] the function 

D,(t) =D q (t;a) = te^ = t(l + i(l - q)at q ~ l )^i , (29) 

where a is a fixed real number. Obviously D q (t) is continuous and differentiable in the interval 
(0, 1). In accordance with the classical Lagrange average theorem for any tx,t2, < t\ < ti < 1 
there exists a number £*, t\ < < £2 such that 

L» ? (ti) - D q (t 2 ) = D' q (U)(ti - t 2 ), (30) 

where D q means the derivative of D q (t) with respect to t. 

Consider the following Cauchy problems for the two Bernoulli equations 

/ 1 ia(l — q) „ , . . . 

V - ~ t V = 1 t V , 3/(0) = 0, (g > 1) (31) 

and 

y' = ^y 9 , i/(0) = (io)^5. (g<i) (32) 

It is not hard to verify that D q (t) is a solution to the problems (|3ip and (|32p with a = x£ for g > 1 
and q < 1 respectively. 

Lemma 3.3 For DJt), g 7^ 1, i/te estimate 

\D q (t;a)\ < C(l + MT^ 1 , * G (0,1], a G #\ (33) 
holds, where C does not depend on t. 
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Proof. Let q < 1. Then it follows from ([32]) that 

\ y '(t)\ = rv| = I4 at? "> = (e^ 2 ' 2 ^ 15 )* < 

C(l + |o|)^5, i E (0,1]. 
Now assume q > 1. Then ()29[) and (|3ip imply 

< t" 1 !!/ + ia(l - g)2/ 9 | = |ef _1 + ia(l - g)(ef ">| = 

|1 +io(l - g)* 9 " 1 !"^ 1 < C(l + |a|) _ ^T, t€ (0,1]. 

Theorem 3.4 Let 1 < q < 1 + oz^n if d >2 and q > 1 if d = 1. Let a sequence of random vectors 
Xjy be weakly q-convergent to a random vector X. Then Xn is q-convergent to X . 

Proof. Assume Xn with associated densities fx is weakly g-convergent to X with an associated 
density /. The difference .F ? [/jv](£) — -^g[/jv](0 can be written 

F q [fN](0-F< I [fN](0= I (D q (f N (x))-D q (f(x)))dx, (34) 
where D q (t) = D q (t;a) is identified in ([29]) with a = x£. It follows from ([30]) and ([33]) that 

i^[/jv](o - -Wjvkoi < c / i(i + wr^ Un(x) - f( x )) \dx, 

JR d 

which yields T q [f N ](0 - .F, [/*](£) for all £ G 

Theorem 3.5 Xei 1 < q < 1 + (^TTJ if d > 2 and 1 < g < 2 i/rf = 1. Assume a sequence of 
random vectors Xn with the associated densities /at is q-convergent to a random vector X with the 
associated density f and Fg [/](£) is continuous at £ = 0. Then Xn is weakly q-convergent to X. 

Proof. Now assume that /at tends to / in the sense of g-convergence. It follows from Lemma 
13.21 that the corresponding sequence of induced probability measures fi^ = PiX^ 1 ) is tight. It 
is well known [151 [T6] that if the sequence of associated probability measures of Xn is tight then 
X]\[ =>■ X and the density of X coincides with /. Lemma [XT] yields Xjj X. 

4 MULTIVARIATE g-CENTRAL LIMIT THEOREMS 

4.0.1 Direct multivariate q-CLT 

A sequence of random vectors X^,N = 1,2,..., is said to be (q d L 1 — q) -independent if, for every 
N=2,3,..., 

r qd jxi + ... + xm) = ® q - ®,^,ra(o, (35) 

where = and X£ = X k - n q (X k ), k = 1,2, .... 

The relation ([35]) can be rewritten as follows. Denote by /at the density of X\ + ... + X C N . Then 
(q'Li — g)-correlation means 

/ e¥ <8u Ma:)rfx = ^[/x f ](0® g .--® a ^[/xd(0. ^ = 2,3,.... (36) 

For g = 1 the condition (|36p turns into the usual independence of identically distributed ran- 
dom vectors X±, ...,Xat- If g 7^ 1, then (gl 1 — g)-independence describes a special type of global 
correlation. 



S 



Theorem 4.1 Assume 1/2 < q < 1 + 2/d. Let Xi, X^f , ... be a sequence of {q d _ l — q) -independent 
and identically distributed d-dimensional random vectors with a finite q-mean \i q = (fi q (Xi), 
and a finite (2q — l)-covariance matrix ^2q-i with the entries cr q ^j = n 2q —x[(Xi — ^q,i){^j — A*gj)]- 
Then Zjq = Xi+ -- + x n i s q d ^ 1 - convergent as N — > oo to a d-dimensional q d _ 1 -normal 

(u 2q -iN) 2 ( 1 + d - d i) 



distribution with the covariance matrix /?X 2(? -i where 







9_1 



4q 

V / 



,q = z d (ql 1 ). (37) 



Proof. Noting that ^ q {X\ — fj, q (Xi)) = 0, one can assume that X\ is centered, that is \x q = 
(0, 0). Let / be the density associated with X\. Using the asymptotic expansion of ^-exponential, 
namely e a q = 1 + a + |o 2 + o(a 2 ), x — ► 0, we obtain 

= L m (i+^oimr 1 - q -(x,if[f(x)f^ +O ((x,o 2 [/(x)] 2 ^- 1 ))) dx = 

l + iv q {i^ q )-{q/2)v 2q ^ 2q ^i) + o{\i\ 2 l |e|-0, (38) 
where (x, £) = sci£i + ••• + and 

= J2^Vq,j = and (E 2 ,-if,0 = I] °"2 9 -lj,feCjCfc- 
j=l i,fe=l 

Hence, for -F 9 [/](£), we have the estimate 

WKO = i - (g/2)^ ? - 1 (s 2? _ 1 e,o + (|e| 2 ), lei - o. 

Further, it is readily seen that, for a given random variable X and a real p > 0, the equality 
F q [pX}(0 = T q [X](p 1 + d - d ^) holds. This yields T q ((u 2q ^NyWT^) X{) = jrj/](^=L=). 

Moreover, it follows from the {q~i — (^-independence of Xi,X 2 ,... and the associativity of the 
g-product that 

F g _ 1 [Z N }(0=F q [f}( ,J )®g-<V 7 g[/]( ,J ) (iV factors). (39) 
Hence, making use of properties of the g-logarithm, from (|39p we obtain 

-|(E 2 ,_,{,{)+o(l), JV-oo, (40) 

locally uniformly by £. 

Consequently, locally uniformly by £, 

lim .F (Z N ) = e -(«M^-^0 e g (41) 

N— >oo y -l 

Thus, Zjv is ^-convergent to the random vector Z, whose (/^-Fourier transform is e q w 2 )^ 2 '- 1 ^). 



9 



In accordance with Corollary 12.31 Z is g^-normal, and hence, its density is a ^Gaussian 
with some covariance matrix S and (3, i.e. 

r q *jz}(0 = r qd jG qd jn- x )\ = e -w^-^\ 

Let us now find £ and (3. It follows from Corollary 12.31 (see (|13p ) that £ = £<•>„- 1 anc ^ ( see CSD ) 

2 ~ 8/3 
Solving this equation with respect to j3 we obtain (f37|l . 



Theorem 4.2 Assume a sequence qf,k = 0, dbl, ... is given as (T7\ ). Let Xjf, N = 1,2,..., be 

a sequence of (qf_ 1 — qf) -independent for some k G Z and identically distributed d-dimensional 
random vectors with a finite qf-mean and a finite second (2qf — l)-variance T, 2q d_ l . 



ThenZ N = D N (qZ)(X 1 + ... + X N - Na d ) with D N (q%) = (AS d ~i) 2(2 ~^ is Qk-r convergent 
to a qf^^normal distribution with the covariance matrix 0k^2q d -\ w ^ ere 



([2-d(qU-l)]C 2 ^~ l) Y d -< 

\ ) 



(42) 



Remark 4.3 The formulation of Theorem \4-S\ changes if we change the definition of q -independence. 
The current formulation, where the definition of {qt-i — qk) -independence is given by \35\). corre- 
sponds to the first raw of Tabled The second and third raws reflect the cases when q-independence 
is replaced by (q^ — q^)- and (qk-i — <Zfc-i) -independence respectively. Consequently, the outcome, 
that is, the type of convergence and the parameter (3k of the attractor are also changed. In all cases 
the attractor is q'^_ 1 - Gaussian (see $5$) with different values of parameter (3 = (3^. Note that all 
three cases generalise the classic multivariate central limit theorem, recovering it if q^ = 1. We note 
also that, in the case d = 1, the outcome of the third raw of Tabled recovers the result presented in 



4.0.2 Sequential multivariate q-CLT 

In this section we consider random vectors X^ with symmetric densities. By definition a sequence 
of random vectors Xn, N = 1, 2, is said to be (Qfc)-independent if 

^oJXx + ... + X N ](0 = ^qJXl](0 ... (8,^ FQ k [X N ](S),VN = 2,3, (43) 

where Tq h is d-dimesional Qfc-sequential Fourier transform and Qk = {qk, qk+i, Qk+d-l}- 
The relation (|43f) can be rewritten as 

^QM(0=^fxA(0^ k+d ---® qk+d ^Qdfx N ](0, N = 2,3,.... (44) 

through the joint density function /jv of X\ + ... + Xjy and density functions fx 3 of Xj. Let f{x) 
be a density. Identify consequently the functions 30 -,x d ), g\{x\, ...,x d -i,U), ■■■ ,9d(£,i, — 

9o( x ) = f( x ), 
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Table 1: Interrelation between the type of correlation, convergence and the parameter of the 
attractor in the multivariate q-CLT. 





Type of correlation 


Convergence 


Gaussian parameter 


1 


F qi jX + Y] = F qi [X] ® qt F qt [Y] 


qf_i — convergent 


Pk — 


q k-l 


\ 1+d - d <-i 

) 


1 4,- 


2 


F q$ [X + Y] = F qt [X] ® qi F qi [Y] 


qf. — convergent 


(3 k = 




i 


4q d 

I 


3 


F qt jX + Y}=F gi jX] ® qt F gi jY] 


qf_i — convergent 


Pk = 


q k-l 


i 

1 + d _ dq d_ i 





gi(xi,...,Xd-i,tid) = 




gj(x±, ...,Xd-j,td-j+ii —,Cd) 



f(*Klf£f k+d - L ' ■■.•■^r~ 3+1W ~ li dx d „ j+1 ...dx d ,j = 2,...,d. (45) 



lR3 



It is not hard to see that 



m<f£ 9 f k+d ~" 1 ■ - ■ ^ ll9d - l]qk ~ 1 dx 1 ...dx d . (46) 



Denote by F^-{x\, ...,Xd-j) the marginal density of f(x), i.e. 



F^Jxi, ...,x d _j) = / f(x)dxd-j+i...dx d 

JR3 



It folows from (|45|) that 

gj(xi, ...,Xd-j,£ d -j+i, .■.,&) = 

Ff^xu ...,x d _j) + o(|&- i+ i| + ... + |&|), + + \£d\ -» 0. (47) 

Further denote 

, x^,_J{x) 

lR d 



o) = f x d _ J f(x)f[[F^ l (x 1 ,...,x d „ l )r+^- l dx, j = l,...,d. 

JR d J 
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Theorem 4.4 Assume 1/2 < q < 1 + 2/ d. Let X±, Xjy, ... be a sequence of Qk-independent and 
identically distributed d- dimensional random vectors with a symmetric density f(x) and a finite 
a], j = l,...,d. 



ThenZ N = (a N ,Y N ), whereY N = Xi+...+X N and on = (a\, aa), aj = (^q k+ j_iNaj) 2 -ik+j-i } 
is q^-convergent as N — > oo to a d-dimensional q^-i-normal distribution with a diagonal covariance 
matrix. 

Proof. It follows from (06]) and flUD that 

(0 = 1 " ^Eft+i-i^J + o(\Z\ 2 )- 

3=1 

Further, it is readily seen that, for a given random variable X and a real d-dimensional vector p > 0, 
with positive coordinates pj, j = l,...,d, there holds FQ k [(p, X)](£) = ^Q k [ X ](Pi 9fc+d_1 £i> p 2 d ~ qk Cd)- 
It follows from this relation that Tq k [(ary, Xi)] = •7 r Q i .[/](;^)- Moreover, it follows from the Qk- 
independence of X\,Xi, ■■■ and the associativity of the ^-product that 

Fq^nKO = ^QJ/K^) - ®« k+d ^QJ/K^) (^factors). (48) 



Hence, making use of properties of the q-logarithm, from (|48H we obtain 

\n qk+d T Qk [Z N U) = N \n qk+d T Qk [f]{^=) = N\n qk+d {l - ^ ^±^M + (^)) = 

-£^^^ + (l),iV-oc, (49) 
i=i 

locally uniformly by £. 

Consequently, locally uniformly by £, 

J im fQ k ( z N) = e qk+d J=1 2 3 . (50) 

TV— >oo 

Thus, ^]\r is Qfc-convergent to the random vector whose Q^-Fourier transform is e gfc+d J=1 ~ J . 
Now we show that for Z there exists ^fc-Gaussian with some covariance matrix £ such that 

*O_*[Z](0 = ^Q fc [G 9fe (E;x)] = e,^ 1 2 

In accordance with Corollary 12.91 the covariance matrix E is diagonal with the diagonal entries 
f3j, j = l,...,d„ uniquely determined from the system of equations 



4a d _ i /3 j 2 

where 7,- and a,-, j = 1, d, are given by (|25|) . The proof of the theorem is complete. 

Remark 4.5 T/ie sequential version of the q-CLT formulated above is valid for the general case of 
non-diagonal covariance matrices. 
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5 Final remarks 



The q-triplet and scaling rate. 

In nonextensive statistical mechanics, the g-triplet (q S en> Qreh Qstat) has been introduced [T7], where 
sen, rel and stat stand for sensitivity to the initial conditions, relaxation, and stationary state, 
respectively. This triplet describes important features of some complex systems, such as the fluc- 
tuating magnetic field of the plasma within the solar wind, as observed in the data sent by the 
Voyager 1 [18]. In paper pQ we considered a different triplet, namely (qk-i, Qk, Qk+i), consisting of 
three succesive members of the sequence {qk, k = 0, ±1, ±2, ...}, corresponding respectively to the 
index of the g-Gaussian attractor (qk-\), the type of correlation, or q-independence, (qk), and the 
scaling rate (qk+i). It appears in fact that there is not just one or other triplet which is relevant, 
but the entire infinite countable family {qk}, each member corresponding to some physical quan- 
tity. The one-to-one association of each member of the family to a concrete physical quantity is 
still an open question; some knowledge is however available (see footnote in [H]). The situation 
seems to follow the same pattern for multivariate processes as well. The direct version of the mul- 
tivariate (/-central limit theorem (Theorem 4.2) provides information about the attractor, the type 
of global correlation, and the scaling rate. More precisely, if the correlation is identified by the 
index qf, then the attractor is a g^^-Gaussian, while the scaling rate is 1/(1 + d— dqf^^ (see, 
the power of N in the scaling factor in Theorem 4.2, i.e., Eq. (|42p . It follows from Lemma 2.5 
that 1/(1 + d — dqf^^ = gf j. Thus, in the d-dimensional case, we have the triplet (<?;?_!,<?;?, Qk+l)- 
If d = 1 we recover the result of [I]. We note that q% is a subsequence of qk = q\ , choosing the 
members numbered nd, n = 0, ±1, and omitting (d — 1) members for each value of n. 

Scale invariance. 

One important concept in statistical mechanics is the concept of exact or asymptotic scale-invariance 
(or scale- freedom). In the present context, we may define this as the case where the joint distri- 
bution (depending on Nd real variables) of a system made of N particles has (either exactly, 
or asymptotically for large N), as marginal distribution, the joint distribution (depending on 
(N — l)d real variables) of a system made of (N — 1) particles (see [HI [H]). We address here 
the study of a similar problem, but related to the dimension d. More precisely, we focus on 
the scale invariance problem for a rescaled process of X\ + ... + where X\, ...,Xn are a se- 
quence of identically distributed Qfe-independent (or (qk-i — (^-independent) random vectors, 

which satisfy the conditions of Theorem 14.41 (Theorem 14.21 respectively). Assume qk > 1 and 

i 

H N (xi,x 2 , ...,x d ) is the density of Z N = (a N ,Y N ), a N = (oi, ...,a d ), aj = (^/q k+ j-iNaj) 2 " 9 *+i-i , 

i 

(or Z N = DN^qOiX! + ... + X N - N^) with D N (q() = (AS^-i) 2(2 ' q " ) )- Then, in accor- 
dance with Theorem 14.41 (Theorem I4T21) . Hn(xi, Xd) — ► Gd, qk _ 1 {x\, ...,Xd) in the sense of weak 
convergence. It follows from Lemma 12.81 and Corollary 12.91 that 

/OO /"OO 
H N (x 1 ,x 2 , ...,x d )dx d = / G d , qk _ 1 (xi,x 2 , ...,x d )dx d = G d -\, qk {xi, ...,x d -\), 
-oo J —oo 

where Gd-i,q k {x\, ...,Xd-\) is a (d — 1) -dimensional ^-Gaussian. Hence, for large N the proba- 
bility density function of the system and its (d — l)-dimensional marginal density have the 
same form of g-generalized Gaussians with the indices qk-i and qk respectively. In other words 
Hn(x\,X2, ■■■,Xd) is, in this specific sense, asymptotically scale invariant. This property enlightens 
the applicability of the present concepts and theorems to many natural and artificial systems, where 
g-Gaussians are observed. 
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